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Abstract 

Context. Rapid rotation modifies the structure of the frequency spectrum of pulsating stars, thus making mode identification 
difficult. 

Aims. We look for new forms of organisation for the frequency spectrum that can provide a basis for mode identification at 
high rotation rates. 

Methods. Acoustic modes in uniformly rotating polytropic models of stars are computed using a numerical code that fully takes 
the effects of rotation (centrifugal distortion and Coriolis acceleration) into account. All low-degree modes, ^ = to 3, with 
radial orders n = 1 — 10 and 21 — 25 for A'' = 3 polytropic models and n = 1 — 10 for A'^ = 1.5 polytropic models are followed 
from a zero rotation rate up to 59 % of the break-up velocity. 

Results. We find an empirical formula that gives a good description of the high-frequency range of the computed acoustic 
spectrum for high rotation rates. Differences between this formula and complete eigenmode calculations are shown to be 
substantially smaller than those obtained with a S'^'^ order perturbative method valid at low rotation rates. 
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1. Introduction 

Asteroseismology has provided a way of probing stellar 
interiors based on the interpretation of observable stellar 
pulsations. In order for such interpretations to be suc- 
cessful, it is necessary to identify pulsation frequencies by 
correctly associating them with theoretically calculated 
pulsation modes. This is important because the geometry 
of a pulsation mode determines which regions of the star 
it probes. One way of doing this is by using regular pat- 
terns which occur in stellar frequency spectra. While being 
quite successful in the case of the Sun and a number of 
slowly rotating stars (Michel, 2006), mode identification 
based on pattern recognition has proved to be very dif- 
ficult in rapidly rotating stars {e.g. Goupil et al., 2005). 
One of the basic reasons for this is that a proper under- 
standing of the structure of the frequency spectrum has 
not yet been achieved in the case of such stars. Up to now, 
mainly perturbative methods, valid at low rotation rates, 
have been used to evaluate the effects of rotation on pulsa- 
tion modes and their frequencies. While providing a useful 
context in which to interpret pulsations of slowly rotating 
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stars, they are unable to correctly predict the structure of 
the frequency spectrum in rapidly rotating stars. This is 
particularly clear in Fig. 5 of Reese et al. (2006), in which 
perturbative calculations of frequencies are compared with 
complete eigenmode calculations for a polytropic stellar 
model rotating at 59 % of the critical angular velocity. As 
a result, many stars remain out of reach for asteroseismol- 
ogy. This mainly concerns early-type stars such as S Scuti, 
which can reach projected equatorial velocities {v ■ s'mi) 
of 200 km.s^^ (Rodriguez et al., 2000), /3 Cephei stars 
{v ■ sini < 300 km.s"\ Stankov & Handler, 2005), and C 
Oph stars (v ■ sin j < 400 km.s^^, Balona & Dziembowski, 
1999). Interestingly, the star ( Oph (w-sini = 380 km.s~^) 
has been observed by MOST and WIRE, thus revealing 
the presence of 19 different pulsation modes (Walker et al., 
2005, Bruntt, private communications). 

Recently, Lignieres et al. (2006a) and Reese et al. 
(2006) came up with a numerical code which overcomes 
the limitations of perturbative methods and enables one 
to accurately evaluate the effects of rapid rotation on stel- 
lar pulsations. By analysing their results, it is possible 
to gain a better understanding of the structure of the 
frequency spectrum at rapid rotation rates. In the case 
of spherically symmetric stars, the origin of the regular 
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frequency patterns can be attributed to the integrabil- 
ity of the ray dynamics which asymptotically describes 
the acoustic wave propagation (Gough, 1993). However, 
in the case of rapidly rotating stars, acoustic ray dynam- 
ics is no longer integrable (Lignicres et al., 2006b), and it 
is an open question whether or not patterns will appear in 
the spectrum of frequencies. Nonetheless, Lignieres et al. 
(2006a) have shown that for axisymmetric modes calcu- 
lated without the Coriolis force, it is possible to obtain 
asymptotic patterns in the spectrum of frequencies. 

In this paper, using again stellar polytropic models and 
a complete non-perturbative computation of their axisym- 
metric and non-axisymmetric acoustic modes, we show 
that regular frequency patterns are present in the acoustic 
spectrum and that they can be used to identify modes. 

2. An empirical formula 

In order to obtain frequency patterns for rapidly rotat- 
ing stars, we start off with the set of pulsation frequencies 
computed in Reese et al. (2006) plus some additional fre- 
quencies. The n = 1 - 10, 21 - 25, £ = - 3, m = to £ 
modes for N = 3 polytropic models, and the n = 1 — 10, 
£ = — 3, m = —itoi modes for N = 1.5 polytropic mod- 
els were followed from a zero rotation rate up to 59 % of 
the break-up velocity. This enabled us to label the modes 
in rapidly rotating stars based on their correspondence 
with modes in non-rotating stars. We looked for patterns 
in this frequency spectrum and found that the frequencies 
approximately obey the following empirical formula in a 
corotating frame: 

i^n,e,m = ^nn + Ai£ + Am\m\ + a"^ (1) 

In order to obtain frequencies in an inertial frame, one 
needs to add the geometrical term —mQ. 

This formula is a generalisation of Eq. (42) of Lignieres 
et al. (2006a) by including the Coriolis effect and the 
case of non-axisymmctric modes. The terms A„, Af, Am 
and a^, which depend on the stellar structure, have been 
computed here for different rotation rates. The param- 
eter a"*" (resp. a~) is an additive constant for symmet- 
ric, £ + m even (resp. antisymmetric, £ + m odd), modes 
with respect to the equator. We note that the term for 
non-axisymmetric modes depends on \m,\ rather than to. 
Indeed, at rapid rotation rates or high radial orders, the 
effects of the Coriolis force become negligible in compari- 
son with the effects of the centrifugal force, as can be seen 
in Fig. 6 of Reese et al. (2006). The centrifugal force leads 
to frequency shifts which do not depend on the sign of m. 
In order to know how these shifts scale with |to|, we plot 
in Fig. 1 the differences {oJn,e,m - '^n,e,m-2.sg(m)) /2 as a 
function of the radial order. We increment the azimuthal 
order by 2 rather than 1 to insure that the two modes 
have the same parity. As shown in Fig. 1, these differ- 
ences converge towards a limit which slightly depends on 
|to| as the radial order increases. This justifies using a lin- 
ear dependence on \m\ as a first approximation. The fact 
that these increments only converge at high radial orders 



shows that Eq. (1) describes an asymptotic behaviour of 
the computed acoustic modes. 
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Figure 1. Frequency differences A™ = "' — """^ as 
a function of the radial order n for the A'^ = 3 polytropic mod- 
els. The frequencies have been scaled by f2,. For high radial 
orders, these differences converge to a limit which does not 
depend on ^ or m. 

In Table 1, we give the values of A„, A^, A^, a+ 
and a~, scaled by Cl^,, for N = 3 polytropic models. 
These are given for different rotation rates scaled either by 

= ^GM/Rl^^ or by VLk = ^GM/R^^, the Keplerian 

break-up velocity, where i?eq (resp. i?poi) is the equatorial 
(resp. polar) radius. These values were obtained by calcu- 
lating various frequency separations and averaging them 
for radial orders 21 < n < 25. At zero rotation, they agree 
pretty well with the theoretical values, A„/n* = 1.238 
and Af/^li, = 0.619 corresponding to the asymptotic de- 
scription of low degree and high order acoustic modes in 
non- rotating stars (Mullan & Ulrich, 1988). The last line 
contains the same parameters but calculated for frequen- 
cies in which the Coriolis force has been neglected. The 
fact that it is essentially the same as the line before shows 
that the Coriolis force plays almost no role in Eq. (1) nor 
ultimately in the frequency spectrum for sufficiently high 
radial orders. 
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0.111 


0.110 


1.239 


0.515 


-0.031 


1.770 


1.843 


0.186 


0.181 


1.234 


0.395 


-0.031 


1.727 


1.942 


0.262 


0.249 


1.229 


0.279 


-0.036 


1.680 


2.006 


0.339 


0.311 


1.220 


0.191 


-0.064 


1.614 


2.020 


0.419 


0.368 


1.194 


0.189 


-0.142 


1.666 


2.068 


0.502 


0.419 


1.154 


0.201 


-0.182 


1.661 


2.036 


0.589 


0.461 


1.102 


0.194 


-0.185 


1.598 


1.953 


0.589 


0.461 


1.102 


0.194 


-0.185 


1.587 


1.943 



Table 1. Coefficients for Eq. (1), for A'^ = 3 polytropes. 



An important consequence of the values given in 
Table 1 concerns the small frequency separation. In non- 
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rotating stars, the so-called small frequency separation 

u>n+i.e.m ^ ^n,e+2.m gocs to zcro In the high frequency 
limit because the ratio A„/A^ is 2. For rapidly rotating 
stars, Table 1 clearly shows that this ratio departs from 2 
by taking on larger values. This, of course, invalidates the 
use of the small frequency separation as a mode identifi- 
cation scheme. 

In Fig. 2, we plot A„ scaled with different quantities 
as a function of As can be seen from the figure, A„ is 
roughly proportional to sjGMjV where M is the mass 
and V the volume of the star. This suggests that A„ may 
be related in a simple way to the mean density of the star. 
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Figure 2. The large frequency separation A„ as a function of 
the rotation rate fi. A„ has been scaled with different quan- 
tities so as to show which other stellar quantity best matches 
this frequency separation. As can be seen, A„ is roughly pro- 
portional to the square root of the mean density of the star. 



3. Accuracy of the formula 

We now turn our attention to evaluating the accuracy of 
the formula (1). In Fig. 3, we plot the correspondence 
between a spectrum based on Eq. (1) and the frequencies 
given by the eigenmode calculations for = 0.59 fix. As 
can be seen, the two sets of frequencies match pretty well, 
especially for high radial orders. This represents a drastic 
improvement over 3'^'^ order perturbative calculations (see 
Fig. 5 of Reese et al. 2006). 

1 /2 

Table 2 contains the mean quadratic error {ey^ both 
for 3'''^ order perturbative calculations and empirical ones. 
These errors as well as the probability p of mode inversion 
(explained below) have been calculated for modes with 
21 < n < 25 except for the last line where 1 < n < 10. 
Physically, an error of 0.0036 A„ at f2 = 0.59 fix corre- 
sponds to 0.16 /iHz for an M = 1.9 Mq, i?poi = 2.3i?0, 
N = Z polytropic star. The errors on the last line can be 
made smaller by calculating the increments A„, A^ etc. 
using radial orders 1 < n < 10, which in practise observers 
are likely to do. 

Another useful quantity to appreciate the reliability 
of an approximate formula at mode identification is the 
probability "p of inverting the identification of two ran- 
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Figures. Comparison between frequencies based on Eq. (1) 
and complete eigenmode calculations, both of which are given 
in an inertial frame. I and m are indicated on the figure, and 
the radial orders are n = 1 to 10. This figure contains four 
columns subdivided into two, the left part corresponding to 
Eq. (1) and the right part to complete calculations. The oblique 
grey lines in between show the correspondence between the 
two sets of frequencies. The units are the same as in Fig. 5 
of Reese et al. (2006), which allows direct comparison. As can 
be seen from comparing these figures, Eq. (1) gives a irmcli 
better description of the frequency spectrum than a order 
perturbative formula. 

domly selected frequencies when using such a formula. 
This probability is defined as: 

Number of inversions 



Number of pairs of modes 



where an inversion occurs for modes A 



/ approx. 



•) {ix>A - OJb) < 0. P = 



(2) 

and B when 
1 corresponds 



to a spectrum in reverse order, whereas p = means the 
order is identical. In Table 2, pis calculated for both meth- 
ods, based on frequencies in the inertial frame so as to be 
closer to observations. For example, p = 0.0035 means 
that out of 3160 frequency pairs, 11 were inverted. As 
shown, the empirical formula gives a far better idea of the 
order in which modes appear in the frequency spectrum 
at ri = 0.59 Hk than a 3'''^ order perturbative approach. 

As expected, the perturbative method gives best re- 
sults for low rotation rates, while Eq. (1) is more efficient 
above fl = 0.11 Ok, for high radial orders. This is because, 
by construction, Eq. (1) preserves the regularity of the fre- 
quency spectrum whereas perturbative methods do not. 
Also, comparing the last two lines of Table 2 again shows 
that Eq. (1) works better for high frequencies, whereas 
the perturbative approach works better for low frequen- 
cies, where the centrifugal force has a smaller effect. 

4. Discussion 

The preceding sections show that at rapid rotation rates, it 
is possible to describe the computed frequency spectrum 
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Perturbative Asymptotic 





\ / 

A„ 
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\ / 

A„ 
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0.037 


0.0005 


0.0000 


0.0114 


0.0000 


0.111 


0.0370 


0.0054 


0.0141 


0.0044 


0.186 


0.1793 


0.0294 


0.0124 


0.0022 


0.262 


0.4325 


0.0595 


0.0135 


0.0047 


0.339 


0.7436 


0.0826 


0.0174 


0.0041 


0.419 


0.9924 


0.1073 


0.0109 


0.0013 


0.502 


1.0674 


0.1237 


0.0085 


0.0025 


0.589 


0.9885 


0.1497 


0.0036 


0.0035 


0.589 


0.3506 


0.0319 


0.1235 


0.0059 



Table 2. Different measurements of the errors for the pertur- 
bative and empirical methods (for N = 3). 



using a formula similar to the asymptotic one found by 
Tassoul (1980). Quite significantly, table 2 shows that the 
accuracy of this formula increases at high rotations, thus 
suggesting that eigenmodes reach an asymptotic regime 
at lower radial orders for these rapid rotations. 

In the rapid rotation regime, a forthcoming study 
based on ray dynamics (Lignieres & Georgeot, 2008) shows 
that Eq. (1) holds for low-degree, high order modes which 
concentrate at middle latitudes as rotation increases. This 
modification of mode geometry leads to a reorganisation 
of node placement as can be seen in Fig. 3 of Reese (2008), 
characterised by a different set of quantum numbers, (n = 
2n + e,l = tiiistil^ ^ — jj^-j ^here e = {I + m) mod 2. 
Based on these quantum numbers, Eq. (1) then takes on 
the following form: 

uJn,e,Tn = ^nn + Aii+ Am\m\+ a"^ (3) 

where A„ = A„/2, At = 2Ai>, A^ = Ae + Am, a+ = a+ 
and — a^+A^— A„/2. With the numerical values from 
Table 1, we find that q+ ~ . The same also applies for 
the iV = 1.5 polytropic model. This suggests that the true 
asymptotic formula is closer to Eq. (3) in which a"*" = 
a~ = a. Applying this new formula leads to errors which 
are only slightly larger in spite of the fact there is now one 
less free parameter. 

Equation (3) also accentuates the |to| dependence of 
the azimuthal term since A^ is of the same order of mag- 
nitude as the variation on A„. This suggests that the 
azimuthal term A„ expressed in these new quantum num- 
bers should take on a different form, and calls for further 
investigation. But as can be seen from Table 2, Eqs. (1) 
and (3) remain effective at identifying low-degree (observ- 
able) pulsation modes. 

For higher degree modes, Lignieres et al. (2006b) 
pointed out the existence of other types of modes, namely 
chaotic and whispering gallery modes, which correspond, 
respectively, to intermediate and high values of i at zero 
rotation. As shown in Lignieres & Georgeot (2008), the 
frequency spectra of these modes have different organisa- 
tions. 



Finally, it must be underlined that the results pre- 
sented here are based on polytropic stellar models. 
Naturally, the question arises whether or not these pat- 
terns still exist in the case of more realistic models. In the 
non-rotating case, sharp chemical gradients lead to the 
break-down of the assumptions behind asymptotic analy- 
sis. Nonetheless, rather than removing the equidistant fre- 
quency pattern, these lead to a periodic component which 
is added to the asymptotic formula {e.g. Gough, 1990). If 
the rapidly rotating case behaves similarly, then it may be 
expected that these sharp chemical gradients also merely 
perturb the asymptotic frequency pattern without remov- 
ing it altogether. 

5. Conclusion 

The formulae (1) and (3) that we presented are the sig- 
nature of an asymptotic regime of high order acoustic os- 
cillations in rapidly rotating polytropic stars. As opposed 
to perturbative methods, which are valid at low rotation 
rates, they give an accurate description of the frequency 
spectrum by preserving its basic structure and conse- 
quently provide a basis for mode identification schemes. 
This is a timely result as seismology space missions are col- 
lecting large data sets of unprecedented quality on rapidly 
rotating stars. 

A key issue for future theoretical studies will be to re- 
late the seismic observables A„, A^ and Am to the physical 
properties of the star. Acoustic ray dynamics combined 
with semi-classical quantization methods is expected to 
play a crucial role in this context. 
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